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ABSTRACT 


The pulsar emission mechanism in the gamma-ray energy band is poorly un¬ 
derstood. Currently, there are several models under discussion in the pulsar 
community. These models can be constrained by studying the collective prop¬ 
erties of a sample of pulsars, which became possible with the large sample of 
gamma-ray pulsars discovered by the Fermi Large Area Telescope (Fermi-LAT). 
In this paper we develop a new experimental multi-wavelength technique to de¬ 
termine the beaming factor (/q) dependance on spin-down luminosity of a set of 
GeV pulsars. This technique requires three input parameters: pulsar spin-down 
luminosity, pulsar phase-averaged GeV flux and TeV or X-ray flux from the as¬ 
sociated Pulsar Wind Nebula (PWN). The analysis presented in this paper uses 
the PWN TeV flux measurements to study the correlation between Jq and E. 
The measured correlation has some features that favor the Outer Gap model over 
the Polar Gap, Slot Gap and One Pole Gaustic models for pulsar emission in the 
energy range of 0.1 to 100 GeV, but one must keep in mind that these simulated 
models failed to explain many of the most important pulsar population charac¬ 
teristics. A tight correlation between the pulsar GeV emission and PWN TeV 
emission was also observed, which suggests the possibility of a linear relationship 
between the two emission mechanisms. In this paper we also discuss a possible 
mechanism to explain this correlation. 

Subject headings: gamma rays: genera, pulsars: general, radiation mechanisms: 
non-thermal, stars: neutron. 
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1. Introduction 

Currently, there are several GeV pulsar models being discussed in the pulsar commu¬ 
nity. These models can be constrained either by studying individual pulsars in detail or 
alternatively using the collective properties of a sample of pulsars. The large sample of 
pulsars discovered by Fermi Large Area Telescope (Fermi-LAT) provides a good place to 
study the collective properties of GeV pulsars. 


The GeV luminosity as a function of other pulsar parameters is a fundamental quantity 
which models must predict. However, the potential utility of pulsar luminosity is limited by 
two factors: inability to measure the beaming factor (/n) (also called the beam correlation 


factor) (Watters et al. 2009) and imprecise distance measurements. The factor /q provides 
the correction to extrapolate the observed phase-averaged flux from the earth line-of-sight 
to the full sky flux for a given beam shape. It is an essential factor needed to convert 
observed fluxes to luminosity: 

Lp = ATid^fnF, ( 1 ) 


where Lp is the luminosity, d is the distance to the pulsar from earth and F is the phase 
averaged flux measured at earth. Since fn is a model dependent parameter, luminosity 
calculations are also model dependent. Therefore, one option to constrain GeV pulsar 
emission models is to use the collective properties of the luminosity distribution but the 
uncertainty on distance measurements degrades the accuracy of the luminosity distribution 
This issue can be resolved by studying the ratio of the flux from pulsars to that of their 
associated Pulsar Wind Nebulae (PWNe), which is a distance independent parameter. 


This paper uses the ratio between a pulsar’s GeV flux and the TeV flux from its 
associated PWN as the first application of this method. Using this ratio we obtain the 
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dependence of /q on spin-down luminosity E for a sample of pulsars. This allows us 
to compare the experimentally measured dependence of /q on E with the theoretical 


expectation of four 7 -ray pulsar gap models (Pierbattista et ah 2012). 


2. The sample of pulsars and their associated PWNe 


Recently the Fermi Large Area Telescope (Fermi-LAT) produced its second pulsar 
catalog ( Abdo et al.||2013 ) with 117 high-confidence 7 -ray ( > 0.1 GeV) pulsars. In addition, 
TeV 7 -ray observatories, such as Milagro, VERITAS and H.E.S.S., have measured TeV fluxes 
coming from PWNe. From a literature survey we found fourteen GeV pulsars in the Fermi- 
LAT pulsar catalog for which the associated PWNe were also measured by TeV observatories. 


Table [T] and summarizes the properties of these 14 objects. This analysis uses the 
values of spin-down luminosity (E), distance to the pulsar (d) and phase averaged flux 
in the energy range of 0.1-100 GeV (Gioo) reported in the Fermi-LAT Second Pulsar Gatalog. 


2.1. Discussion of TeV PWNe Measurements 


The integrated energy flux around 35 TeV of the associated PWNe {Etsv) are listed in 
column 3 of Table All Milagro TeV measurements in this column are derived from Table 


1 in Abdo et ah (2009a). Hereafter, we will refer to Abdo et al. (2009a) as the Milagro 


OFGL search. In that publication the Milagro collaboration performed a targeted search for 


galactic sources in the Fermi Bright Source List, which is also known as OFGL (Abdo et 


al. 2009b). The Milagro OFGL search found TeV emission coincided with 14 Fermi bright 
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sources. Among these 14 sources, 9 have spatial associations with pulsars in the Fermi-LAT 
Second Pulsar Catalog. The Milagro OFGL search paper reported the differential photon 
ffux at 35 TeV I 35 Tev) assuming a Spectral Energy Distribution (SED) of E“^'® without 
a cutoff. The authors argue that the ffux calculated at 35 TeV has the least dependence of 
the calculated ffux on the true spectrum. Because Milagro does not report the full SED for 
these pulsars, but only the differential photon ffux at 35 TeV, we calculated the integrated 
energy ffux over a 1 TeV band around 35 TeV with a SED of E“^'® and used the Milagro 
ffux uncertainties. 


The Milagro OFGL search paper mentioned that TeV emission might come from the 
pulsar and/or from the associated PWN. However, it is very unlikely to get a significant 
contribution from the pulsar to the TeV ffux measured by Milagro. The best example of this 
is the Grab pulsar, which is the brightest TeV object measured by Milagro. The VERITAS 


Gollaboration (VERITAS Gollaboration et ah 2011) observed pulsed 7 -rays in the energy 


range of ~ 100 GeV to ~ 200 GeV. The measured energy spectrum is well described by a 
simple power law, without a cut-off, by: 


dN , 

— = i^A.2±0.Qstat±lAsyst 


2.4 


X 10 


-11 


E 


150 GeV 


—3.8di0.5s£ati0.2syst 


TeV^cm-V^ (2) 


An extrapolation of this energy spectrum gives a differential photon ffux of 
4.2 X 10“^° photons TeV“^cm“^s“^ at 35 TeV. This is 0.003% of the TeV ffux ob¬ 
served coincident with the Grab pulsar by Milagro. In addition, a theoretical model 


proposed in Aharonian et ah ( 2012 ) predicts a sharp cut-off below ~500 GeV, so the 


extrapolated ffux at 35 TeV from the pulsar might be even lower. These considerations lead 
us to conclude that the TeV emissions observed coincident with pulsars come predominantly 
from their associated PWNe. 
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We performed a literature search for PWNe measured by Air Cherenkov telescopes 
in the TeV band. We found H.E.S.S. SEDs for Eve other PWNe, which are also listed in 
Table [T] and 1^ Crab, Vela, K3 in Kookabura, MSH 15-52, and G 21.5-0.9. In order to be 
consistent with the Milagro measurements, we integrated the energy flux within a 1 TeV 
energy band around 35 TeV using the H.E.S.S. SEDs. The SEDs of the Crab, Vela and 
MSH 15-52 PWNe were measured in the energy ranges of 440 GeV - 40 TeV, 550 GeV - 65 
TeV and 250 GeV - 40 TeV, respectively. Therefore, their integrated energy flux around 35 
TeV can be obtained without any extrapolation. However, the SEDs of K3 in Kookabura 
and G 21.5-0.9 PWNe were measured in the energy ranges of 200 GeV - 25 TeV and 150 
GeV - 5 TeV, respectively. Therefore, their integrated energy fluxes obtained around 35 TeV 
by extrapolation of their SEDs:ql might be an overestimate, if there is a cutoff below 35 TeV. 


VERITAS has published SEDs of the Boomerang and GTA 1 PWNe. In both cases 
the SEDs of these sources were measured in the energy range of 1-15 TeV. Therefore, as for 
K3 and G 21.5-0.9, the integrated energy flux obtained around 35 TeV by extrapolating 
the SEDs might be an over estimate, if there is a cutoff before 35 TeV. For all H.E.S.S. 
and VERITAS measured PWNe, errors on the integrated flux are estimated by a standard 
Gaussian Monte Garlo propagation of the uncertainties of the SED £t, with the 16*^ 
percentile as the lower error bar and the 84**^ percentile as the upper. 


There are two independent measurements for both Boomerang and the Grab. In each 
case, the measurements agree within experimental errors and differ by less than a factor 
of two. Both measurements for these PWN are shown in the following plots but we use 
their weighted average when doing hts, which does not alter any of the conclusions in this 
analysis. 


Table 1. Properties of a sample of GeV pulsars cataloged in the Fermi-LAT Second Pulsar 

Catalog. 


Pulsar Name 

Association 

logio(^ 

Gioo (10 

Distance 

PSR 


( 10 ^"^ erg s“^)) 

erg cm“^ s“^) 

kpc 


J0007+7303 

CTA 1 

35.65 

40.1+0.4 

1.4+0.3 

J0534+2200 

OFGL J0534.6+220lf 

38.64 

129.3+0.8 

2.0+0.5 

_II_ 

Crab 

_II_ 

_II_ 


J0631+1036 

OFGL J0631.8+1034f 

35.24 

4.7+0.3 

1 . 0 + 0.2 

J0633+1746 

OFGL J0634.0+1745t 

34.52 

423.3+1.2 


J0835-4510 

Vela 

36.84 

906+2 

0.29 + 0.02 

J1420-6048 

K3 in Kookabura 

37.01 

17.0+1.4 

5.6+0.9 

J1509-5850 

MSH 15-52 

35.71 

12.7+0.7 

2.6+0.5 

J1833-1034 

G2L5-0.9 

37.53 

5.9+0.5 

4.7+0.4 

J1907+0602 

OFGL 1907.6+0602t 

36.45 

25.4+0.6 

3.2+0.3 

J1958+2846 

OFGL 1958.1+2848f 

35.53 

9.1+0.4 

< 18.5 

J2021+3651 

OFGL 2020.8+3649^ 

36.53 

49.4+0.8 

io.oli° 

J2021+4026 

OFGL 2021.5+4026f 

35.06 

95.5+0.9 

1.5+0.4 

J2032+4127 

OFGL 2032.2+41221 

35.44 

10 . 6 + 0.6 

3.7+0.6 

J2229+6114 

OFGL 2229.0+61141 

37.35 

25.3+0.4 

0 . 8 +l°i^ 


Boomarang 






Note. — Gioo is the phase averaged flux of the pulsar GeV emission in the 0.1-100 


GeV energy band. 
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TeV flux of the associated PWNe of a sample of GeV pulsars cataloged in the 
Fermi-LAT Second Pulsar Catalog. 


Pulsar Name 

Association 

o 

1 

PSR 


TeV s ^ cm 


J0007+7303 

CTA 1 

1 /I+4.9“ 
-‘-•^-1.1 

J0534+2200 

OFGL J0534.6+220lf 

5.4 ±0.3^ 

_II_ 

Crab 

o q+6.0^ 
^•^-2.0 

J0631+1036 

OFGL J0631.8+1034f 

1.5 ± 0.4^ 

J0633+1746 

OFGL J0634.0+1745t 

1.2 ± 0.4^ 

J0835-4510 

Vela 

16.4+f 

J1420-6048 

K3 in Kookabura 

5.4±{r 

J1509-5850 

MSH 15-52 

6 .2±»f 

J1833-1034 

G21.5-0.9 

0.99±i/ 

J1907+0602 

OFGL 1907.6+0602t 

3.9 ± 0.5^ 

J1958+2846 

OFGL 1958.1+2848^ 

1.1 ± 0.3^ 

J2021+3651 

OFGL 2020.8+3649^ 

3.6 ± 0.3^ 

J2021+4026 

OFGL 2021.5+4026^ 

1.2 ± 0.3^ 

J2032+4127 

OFGL 2032.2+41221 

2.1 ± 0.3^ 

J2229+6114 

OFGL 2229.0+61141 

2.3 ± 0.4^ 


Boomarang 

1 q+15 h 
i.O_i 2 
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Note. — Fn is the integrated energy flux of the PWN 
TeV emission in the 34.5-35.5 TeV energy band. 


[VERITAS Measurement. Energy flux derived by ex¬ 


trapolating the SED.] Reference Alin et al. (2013) 


’[Milagro Measurement] Reference Abdo et al. (2009a) 


’[H.E.S.S. Measurement] Reference Aharonian et al 


(2006a) 


'[H.E.S.S. Measurement] Reference Aharonian et al. 


(2006b) 


'[H.E.S.S. Measurement. Energy flux derived by ex¬ 


trapolating the SED.] Reference Aharonian et al. (2006c) 


^[H.E.S.S. Measurement] Reference Aharonian et al 


(2005) 


^[H.E.S.S. Measurement. Energy flux derived by ex¬ 


trapolating the SED.] Reference|H. E. S. S. Collaboration 


et al. (2007) 


^[VERITAS Measurement. Energy flux derived by ex¬ 


trapolating the SED.] Reference Alin et al. (2013) 


























11 


3. Method 

The ratio of the pulsar GeV luminosity (Tp) to the luminosity of the associated PWN 
(Ltv) can be written in terms of the corresponding pulsar and PWN Flux: 



G^ioo 


fn X 


’AT 


( 3 ) 


Taking their ratio cancels the distance but retains the beaming factor f^, which can be 
written as: 

fn = Fr (4) 

-^n 

where 

Fn /-X 

r = —, (5) 


Gioo’ 


which is the observed flux ratio. This relation is a mathematical identity which is valid for 
each individual pulsar and its associated PWN. This identity however can not be used to 
derive for a given pulsar, because Lp is not measurable without fn. However, we can 
extract information on the E dependence of fn for a selected group of pulsars by using 
models that predict the E dependence of Lp and Ltv. 


Many high energy pulsar models (e.g. Harding (1981); Harding & Muslimov (2002); 
Takata et ah ( 2010| ; Muslimov & Harding (2003) and Harding & Muslimov (2003)) predict 


a power law relationship between Lp and E. 


Lp = kp- E^ (6) 

For a given pulsar, kp is independent of E, but depends on other pulsar properties such as 
the angle between the direction of the magnetic dipole axis and the rotation axis. Both kp 
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and the power q are model-dependent. Later in the paper, we will discuss the implications 
of different choices of q. 


For PWNe 


Mattana et ah (2009) discussed the correlations between E and PWN 


luminosity in TeV and X-ray energy bands. Using H.E.S.S. measurements they showed 
that PWN TeV luminosity is not correlated with E. This observation is consistent with 
the theoretical expectation, TeV photons are generated by the accumulated high-energy 


electrons in PWNe (Mattana et ah 2009). Therefore, for a given ensemble of GeV pulsars 


we can choose a characteristic PWNe TeV luminosity /cjv, independent of E. Mattana et ah 


(2009) also showed that the X-ray luminosity vs E distribution can be htted into a power 
law model. Therefore we can generalize the X-ray luminosity vs E distribution and TeV 
luminosity vs E distribution as; 

LN = kN- E^, (7) 


where m = 0 for TeV luminosity vs E distribution. 


Both of these energy bands are good candidates for applying our method. The work 
presented in this paper uses PWN TeV luminosity. Another analysis with PWN X-ray 
luminosity is in progress. 

When we combine the model expectations in Equation for pulsars and Equation for 
PWNe with Equation]^ we obtain for a specihc pulsar v. 

In log-log space we can rewrite this equation as: 

log/oi = log {j^ + + {q-m) log (Ei). (9) 


In this equation r* and Ei are measurable quantities, but the coefficients kp and k^ are 
unknowns and vary pulsar to pulsar. In this paper we do not intend to measure the of 
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individual pulsars. Instead we intend to obtain the /q dependence on E for an ensemble of 
GeV pulsars using kp and kjq, where kp and kjq are typical values of kpi and krii appropriate 
for our ensemble of pulsars. We make this explicit by defining di as the difference between 
typical values and the pulsar-dependent constants: 

We can rewrite Equation with the parameter di as, 

log/o* - di = log I ^ ) + logTi + {q-m) log (Ei). (11) 

\kN J 

Although di is not measurable for individual pulsars, we can use this expression to obtain 
the dependence of an estimate of (fm) on Ei^ where 


log fm = log fm - di 

= log I I -h log ri + {q- m) log Ei 
kN 


( 12 ) 


Thus di is a correction factor between using typical values and the unknown pulsar-dependent 
values. We will estimate the magnitude of any such effects in section 


This summarizes our method of extracting the E dependence of /q. We now proceed 
to discuss our choices for the constants g, m, kp and kj^j in more detail, and examine how 
well data supports these choices. 


4. TeV PWN Measurements 

First we consider whether pulsar and TeV PWN luminosity exhibit sufficient correlation 
to make it worthwhile to work with their ratio. Figure [T] shows the correlation between 
Lp/fn and Lpf. The blue squares are PWNe measured by Milgro, the red stars by 
H.E.S.S., and the green circles by VERITAS. It appears that the PWNe measurements 
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Fig. 1.— PWN TeV luminosity vs. Pulsar GeV luminosity normalized with respect to 
the beaming factor (/o), dvrd^Gioo- Blue squares are measured by Milagro. Red stars are 
measured by H.E.S.S. Green circles are measured by VERITAS. The two dotted lines have 
a slope of 1 but different arbitrary intercepts. The linear correlation coefficient is R = 0.82. 
The error bars are dominated by the distance measurement uncertainties. 
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of the three experiments are generally consistent. Four of the TeV measurements 
(Ltv = 34.3,35.0,35.2,35.8) have been extrapolated from lower energy, but do not appear 
to be outliers to the general scatter of the distribution. The large error bars on this plot 
are due to the distance uncertainties. Also note that the TeV luminosity (y-axis) spans a 
larger range than the GeV (x-axis), making the error bars appear to be less important for 
TeV luminosities. The dotted lines (drawn to guide the reader’s eye) represent two lines 
of Lat oc ^ with different arbitrary intercepts. Even though the error bars are large a 
reasonable correlation is obtained with a linear correlation coefficient of R = 0.82, which we 
judge to be sufficiently encouraging to proceed. 

Next we examine the correlation between and E. The Ljsi vs E distribution 

for our PWN sample is shown in Figure Again we note that the uncertainty on the 
distance measurements contributes signihcantly to the luminosity error bars and that 
the extrapolated points {E = 35.7,37.0,37.3 and 37.5 ) do not appear to be outliers. 

This distribution has a linear correlation coefficient of 0.09. The small linear correlation 
coefficient suggests that Tat is not correlated with E that concludes PWN TeV luminosity is 
not a function of E. Therefore, one has to expect zero slope for the best £t linear fit for the 
data points. The best fit linear fit for our data points has the slope of 0.03 ± 0.06, which is 
consistent with zero. In summary, we argue that the observations are consistent with the 
theoretical expectation of no E dependence in the PWN TeV luminosity, which we discussed 
in section Therefore, the model value m = 0 in Equation We £t a constant to the 
loglvAr data, yielding 31.6 ± .05, and use this value for the model parameter \ogkp^ = 31.6. 


5. GeV Pulsar Measurements 

Next we return to the power law model for the pulsed GeV emission Lp. We proceed 
with the analysis on the basis of g = because there are several high energy pulsar models 
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Fig. 2.— PWN TeV luminosity vs. spin-down luminosity E of the associated pulsar. Blue 
squares are for PWNe measured by Milagro, red stars by H.E.S.S., and green circles by 
VERITAS. The distance uncertainty contributes signihcantly to the error bars. 
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which predict a power law index near | (e.g. Harding & Muslimov 


Takata et ah 


(2010); Muslimov &: Hardingj (2003) and Harding & Muslimov ( j2003 )), and further, because 
in this paper we compare our results with a sample of pulsars simulated using Lp = kpE^ 
as an underlying model assumption. Later we will discuss the effect of q on the best fit 
parameters in Equation [TOj 


Finally we can select a reasonable value of /cp by using Figure 9 of 
which has an illustrative line for a Lp = kpE^ model with an additional constraint of 
Jq = 1. Taking a suitable point from the line, logF = 39 and \og{Lp/f^i) = 36, we find 
log kp = 16.5. 


Abdo et al. 


(2013 


6. Analysis 


As a summary of Sections and we proceed with our analysis using model 
parameters q — m = \ and logi^ikp/kn) = —15.1. With these model parameters we can 
rewrite Equation [T^ as follows: 

logio hi = logio(ri -Ef)- 15.1 (13) 

The correlation between log^g = yi and \ogiQ Ei = x* is shown in Figure]^ It appears 
that above logF 35, this distribution has a linear correlation with 

y = (-11.04 ± 1.13) + (0.28 ± 0.03) ■ x^, (14) 


where 


yi = log fm 
= log hi - di 

= log ( ^ ) + log ri + {q- m) log Ei 
kN 


Xi = log Ei 


(15) 
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and x^/NDF for this fit is 26.9/14. 

Therefore, we can ht a phenomenological power law model for logE > 35, 


fn = a + b-\ogE, 


(16) 


This power law model is motivated by Pierbattista et ah (2012), and it also appears to be 


in agreement with the data. The y^/NDF (= 26.9/14) for the £t in Equation 14 indicates 
deviations of from the fit line {yi) due to a combination of statistical measurement errors 
and systematical errors. The systematic scatter can be estimated by adding a constant e 
in quadrature with the statistical errors of each yi, such that the chisquared per degree 
of freedom decreased to 1.0. The resulting estimate is e = 0.08, which indicates ~ 20% 
systematic error. The slope of the £t with this added uncertainty is 0.27 ± .06, which is still 
inconsistent with zero. 


The fit residuals represent the difference between the data points and the empirical 
htting model: 

yi-yi = logfm -di-yi (17) 

Thus e characterizes the scale of the differences among the estimated functional form 
(trend) y, and the individual pulsar measurements, effectively considering the systematic 
deviations from the empirical power law for Jq as well as the effects of removing d* from Jq 
(that is using the single pulsar-independent values of kp and /cjv)- We see that although 
these deviations represent the information about specific pulsars compared to the overall 
model (trend), the deviations of individual pulsars from the trend are not so large as to 
invalidate the model extraction of the trend, as this scale (e = 0.08) is notably smaller than 
the variation of y across the range of E. 
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Fig. 3.— Experimentally obtained /q vs. pulsar spin-down luminosity E. Blue squares are 
PWNe measured by Milagro, red stars by H.E.S.S, and green circles by VERITAS. The best 
fit line for the pulsars with E > ergs s“^ is shown as a dotted line. The slope of the 

best £t line is 0.28 ± 0.03. 


Discussion 


7.1. Constraints on y-ray pulsar gap models 


Recently, Pierbattista et ah (2012) studied four gamma-ray pulsar acceleration models. 


They synthesized a pulsar population based on a radio emission model and four y-ray 
pulsar gap models: Outer Gap (OG), Polar Gap (PG), Slot Gap (SG) and One Pole 
Gaustic (OPG). Their model simulations of the correlation between /q and E are shown 
in Figure]^ In all four model predictions, for pulsars with E > 10^® ergs s“^, /n can be 
reasonably ht by a straight line in log — log space and the best hts give the following slopes: 
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mpc = 0.05 ± 0.23, msc = 0.003 ± 0.004, moc = 0.12 ± 0.01 and mopc = 0.026 ± 0.007]^ 
The slopes of the PC and SG models are consistent with zero, and the slope of the OPC 
is very small. The OG model is the only model that predicts a positive slope that is not 
consistent with zero. However, the small number of data points for the PC model have a 
large scatter and while the slope is consistent with zero, the uncertainty on the slope is 
much larger than for the other models. 


We can now compare our data points from Figure with the expectations from the four 
pulsar gap models. First we note some gross characteristics in comparing the simulations 
to our data points. Both the data and the models have log /q near 0, and the scatter of 
individual simulated pulsars about their trend line is not grossly different from the scatter 
in the data, despite the fact that the data extraction used single individual values oi 
and kf), while the simulation used the full information about individual simulated pulsars. 
This similarity is in agreement with the idea that the scatter due to Cj is not so large as to 
lose all information about fn. These simulated data points show a tighter distribution of 
the simulated radio-loud pulsars to the best fit line above log^Q(i7) = 35 while radio-quiet 
pulsars have a wider distribution. Especially in the OG model radio-quiet pulsars deviate to 
low /q values for log^Q(£') < 35. We note that the radio-loud pulsars in our data also have a 
tighter distribution about the best fit line above logiQ(i7) = 35 and the radio-quiet pulsars 


with log]^g(i7) < 35 deviate to low /q values. The authors of Pierbattista et ah (2012) also 
noted that the range of variation of their was less than that of their Lp for the same 
range of E, at least for the SG and OG models. We see a similar trend in our data. 


The experimental /q vs. E distribution has a non-zero slope of 6 = 0.27 ± 0.03 for 
E > 10^® erg s“^. This would tend to disfavor the PC, SG and OPC models, despite the 


^These best fit parameters were provided to us by the authors of 


Pierbattista et ah 


( 2012 ). 
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Fig. 4.— The dependence of the beaming factor /q on spin-down luminosity E for four 


models derived by Pierbattista et ah (2012) a sample of simulated pulsars. This plot was 
reproduced by M. Pierbattista with linear hts for pulsars above E = 10^^ ergs s“^. For 


the original hgure refer to Pierbattista et ah (2012). Red and green markers refer to the 


radio-loud and radio-quiet pulsars, respectively. Black lines refer to the best linear hts for 
the pulsars with E > 10^® ergs s“^. Slopes of the best ht lines are mpc = 0.05 ± 0.23, 
'msG = 0.003 ± 0.004, moo = 0.12 ± 0.01 and mopc = 0.026 ± 0.007. 
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OPC providing the best overall agreement with Fermi-LAT pnlsars among the models 
considered. However, the slope in our data is over twice that expected by the OG model for 
E > 35 ergs s“^. Below that value of E, the expected correlation between /o and E for the 
OG model becomes more dispersed and the /o distribution for a given E has a tail towards 
smaller /q values, especially for radio quiet pulsars. This feature is also consistent with the 
experimentally obtained /n vs. E distribution. The two radio quiet pulsars with E below 
Iqss g-i (^pgj^ J0633+1746 and PSR J2021+4026) have smaller values compared to 
radio loud pulsars. By considering both features above and below E ^ 10^^ ergs s“^, we 
can conclude that our data sample has some features which favor the OG model for pulsar 
emission in the energy range 0.1-100 GeV over the SG and OPG models, though even the 
OG model does not quantitatively match our measurements. However, we cannot reach any 
conclusions about the PG model. 


The discrepancies with the models might be due to the systematic limitations or 


inadequacies of Pierbattista et ah (2012)’s simulations or biases in our data sample. As 


the authors of Pierbattista et ah (2012) mention, results of all four simulated models are 
lacking pulsars visible by Fermi-LAT with E > 3 x 10^^ ergs s“^ and characteristic age 
< 100 kyr, and over-predict the number of low E pulsars. Furthermore they also mention 
that the simulated OG model in particular fails to explain many of the most important 
pulsar population characteristics, including the distributions of period, characteristic age, 
Lp and E. This could certainly affect their Jq vs. E distributions. Therefore, we cannot 
provide tight constraints using this synthetic pulsar population. However, a comparison of 
our results with an improved model simulation could provide tighter constraints. 


A better sensitivity TeV survey such as HAWG promises to help our understanding 
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of PWN visibility selection effects. Our estimates of may also be biased if there is an 
unexpected residual dependence of the PWN TeV luminosity on E hidden in the TeV data 
scatter, or if there is hidden dependence on the selected TeV energy range. The extracted 
slope of the vs E distribution linearly depends on any hidden slope in the PWN TeV 
luminosity vs E distribution. The TeV energy range of a 1 TeV band around 35 TeV was 
chosen because of the Milagro data set. Ideally, we would prefer to do this study with a 
more uniform sample of PWN TeV energy fluxes obtained around the inverse Compton 
peak of each individual PWN. Using X-ray PWN luminosity instead of TeV luminosity 
would also offer complementary visibility selection effects, to allow assessment whether such 
selection effects are important. 


This result also depends on the uncertainties of the theoretically predicted correlation 
between GeV pulsar luminosity and spin-down luminosity, Lp oc E'^, q = ^ + 6. We made 
Figure]^ for Lp oc E^ (5 = 0), because it is natural in several pulsar models to have this 
relation (|Abdo et al. 2010|[) and it is one of the underlying assumptions in Pierbattista et ah 


(2012). However, the slope of log^g(/n) vs. log]^o(-^) distribution depends linearly on 5, as 
seen in Equation 15 Therefore we can write that for E > 10^^ erg s“^ the slope of \ogif^{fQ) 
vs. log]^Q(£') distribution is 0.28 ± O.OO^tat + 5. One can use this expression to determine the 
slope of logiQ^fn) vs. logiQ(E) distribution under different models. For example. 


Takata 


et al. 


(2010) discuss a high-energy emission from the outer gap that expects Lp oc Es. 


Under this model the slope 6 = 0.125 and slope becomes 0.405 ± 0.03stat- Another example 

is 


Muslimov & Harding (2003), which discuss a high-energy pulsar emission model that 


depends on the local magnetic held. In the low magnetic held scenario Muslimov fc Harding 


(2003) expects Lp oc Er,6 = —0.07. Under this model the slope becomes 0.21 ± 0.03stat. 
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7.2. Pulsar GeV emission to PWN TeV emission connection 

The correlation between the pulsar GeV emission and the PWN TeV emission shown in 
Figure [pleads one to suspect a common underlying cause for the two emission mechanisms. 
One property relevant to both emissions is the electron-positron current of the pulsar 
wind {Iwind)- The GeV energy flux from pulsars is thought to be directly related to the 
instantaneous value of Iwind^ because the GeV pulsed emission from the magnetosphere 
is often thought to be produced by curvature emission by the most recently produced 
electron-positron population in the wind. Since the luminosity is roughly proportional to 
the population of electrons and positrons, we can write: 

LpSR GeV OC Iwind- (18) 


Inside PWNe, TeV photons are often thought to be produced by the up-scattering of 
ambient photons by the relativistic electrons and positrons, known as inverse Gompton 
radiation. Therefore, Tat should depend on the relativistic electron-positron population 
and the ambient photon population in the PWN. However, for the relativistic electrons and 
positrons that produce TeV photons by inverse Gompton scattering the typical cooling time 


is larger than the lifetime of pulsars (Mattana et ah 2009). Therefore, the population of 
these electrons and positrons becomes proportional to the integral of Iwind over the pulsar 
lifetime, instead of proportional to the instantaneous value of Iwind- However, we could 
suggest a proportionality between the ambient photon held density (pph) and Iwind- There 
are two different ambient photon helds which could be relevant to the production of TeV 


7 -rays: photons from synchrotron radiation and far-infrared photons (Atoyan & Aharonian 


1996). The density of synchrotron radiation photons in the x-ray energy band is roughly 


proportional to the density of the freshly injected pulsar wind (Mattana et al. 2009). In 
addition far-infrared seed photons can be made by heating the pulsar wind, as described in 


Section 2.2 of Arons (1996). Therefore the pph may be roughly proportional to Iwind- If Pph 
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is proportional to Ljsf, that would yield: 


L]\l OC Pph OC lujind- 


(19) 


Hence, 


L]Sf OC LpsR GeV- 


( 20 ) 


While these considerations are suggestive, a more detailed theoretical study is clearly 
needed to fully understand this correlation. 


8. Conclusion 


We have developed a new multi-wavelength technique to study the collective properties 
of the GeV pulsar beaming factor Jq with respect to the pulsar spin down luminosity E. This 
technique uses a distance independent parameter, to obtain the correlation between 
fn and E. It allowed us to use the pulsars with poorly known distance measurements to 
study Jq. Using this technique we have experimentally obtained the fn vs. E dependence 
for pulsar emission in the 0.1-100 GeV energy band. Under the model assumptions of an 
dependence of GeV pulsed emission but no E dependence of TeV PWN emission, we hnd 
a dependence of /n on E. Our experimentally obtained correlation between /q and E has 
some features which favor the theoretical fn vs. E distribution of the OG model obtained 


by Pierbattista et ah (2012). However, this specihc comparison is limited by the modeling 


uncertainties of Pierbattista et ah (2012)’s simulated pulsar sample. Applying this same 
multi-wavelength method to X-ray data for PWNe is attractive, since it may have a more 
precisely measurable E dependence than the present TeV data. 


Pulsar GeV emission and PWN TeV emission are correlated, with a linear correlation 
coefficient of R = 0.82, although TeV PWN emission has no correlation to E. This observed 
GeV to TeV correlation suggests the possibility of a linear relationship between y-ray 
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emission mechanisms in pulsars and TeV emission mechanisms in PWNe. However, it is 
not possible to explain this linear relationship using the electron-positron populations of 
curvature radiation inside the magnetosphere and synchrotron radiation in the PWN. An 
alternative possibility is a linear relationship between the ambient photon density {pph) in 
the PWN and the pulsar wind current {Iwind)- A more detailed theoretical study will be 
needed to fully understand this correlation. 

In the near future, TeV experiments under development, such as HAWC, CTA, and 
Lhasso, will have greater sensitivity than Milagro. The observed GeV to TeV luminosity 
correlation makes it likely that these observatories will detect PWNe associated with many 
more of the GeV pulsars Fermi has observed, leading to prospects of a higher-statistics and 
higher precision version of this analysis. 


We would like to give our special thanks to the authors of Pierbattista et ah (2012), 
particularly to Dr. M. Pierbattista for producing new hts for the simulated data and Prof. 
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thank Prof. K. Tollefson and Dr. J. Pretz for useful suggestions, and Dr. G. Sinnis and 
Prof. M. Voit for their encouragement to publish this work. This work was supported by 
the National Science Foundation. 
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